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Multiple  scattering  of  classical  particles  on  a system  of  point  scatterers  distributed 
in  space  in  a fractal  fashion  is  considered  in  the  small-angle  approximation.  The 
asymptotic  regime  of  this  process  is  described  by  a fractional  differential  equation 
generalizing  the  ordinary  diffusion  in  the  angle  space  equation.  The  solution  of 
the  problem  is  presented  both  in  terms  of  stable  distributions  and  in  terms  of 
Pox-functions.  Main  features  of  the  obtained  solutions  are  discussed. 


1 Introduction 


Small-angle  x-ray  and  neutron  scattering  experiments  are  efficient  tools  for  study- 
ing the  structure  of  condensed  matter.  Measurements  of  the  light  or  cosmic  rays 
from  distant  sources  play  the  same  role  in  studying  the  large-scale  structure  of  the 
Universe.  From  mathematical  point  of  view,  the  problems  belong  to  the  class  of 
inverse  problems  solved  in  terms  of  the  multiple  scattering  theory. 

The  ordinary  multiple  scattering  theory  assumes  that  the  random  spatial  distri- 
bution of  scatterers  is  a homogeneous  Poisson  ensemble,  i.e.  different  scatterers  are 
placed  independently  of  each  other.  That  is  just  what  guarantees  the  exponential 
character  of  the  free  path  distribution  and  leads  to  the  Boltzmann  master  equation. 
In  the  region  of  large  depths  and  small  angles  the  equation  takes  the  form  of  the 
diffusion  (in  angle  space)  equation  called  the  Fokker- Planck  equation.  Its  solution 
is  nothing  but  the  two-dimensional  (in  angle  space)  Gauss  distribution  called  also 
the  Fermi  distribution.  However,  above-mentioned  assumption  becomes  invalid  for 
fractal  media  characterized  by  long-range  correlations  of  the  power  type. 

One  approach  to  solution  of  the  problem  is  developed  by  S. Maleyev  1 . He 
has  considered  a medium  with  fractal  pores  distributed  uniformly  in  space.  As 
a result  the  superdiffusion  behaviour  of  the  angular  distribution  has  been  found. 
Another  approach  is  developed  in  our  previous  work  2,  where  the  fractal  medium 
is  considered  as  a set  of  identical  point  scatterers  distributed  in  space  in  fractal 
manner  by  means  of  the  Levy  flight  procedure  (see  for  details  3’4,5).  This  model 
reveals  the  subdiffusion  behaviour  of  the  angle  distribution. 

The  present  work  combines  both  of  these  ideas.  We  consider  here  a set  of 
fractal  clusters  distributed  in  space  in  a fractal  manner.  As  a result,  we  obtain  a 
two-parametrical  family  of  angular  distributions.  The  family  includes  the  Fermi 
distribution,  the  subdiffusion  and  superdiffusion  distributions  as  particular  cases. 
They  obey  fractional  differential  equations  and  are  presented  both  in  terms  of  stable 
distributions  and  in  terms  of  Fox-functions. 
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2 Basic  postulates  of  the  small-angle  scattering  theory 

Let  a particle  move  from  the  origin  of  coordinates  along  the  z-axis.  After  the 
scattering,  the  particle  will  be  characterized  by  the  vector  angle  © of  deviation 
from  the  initial  direction  flo  of  the  motion 

© = — Oo- 

The  following  postulates  form  the  well-known  basis  of  the  small-angle  approxima- 
tion in  the  scattering  theory. 

i)  The  particle  is  on  the  z-axis  during  all  its  motion  even  though  the  (small) 
vector  0 differs  from  zero. 

ii)  Vector  0 stays  constant  between  points  of  scatterings  X\ , X1  + A2,  Xj  -|  A2  4 - 
A3, . . . and  undergoes  random  jumps  ©2,  ©2,  ©3, ...  at  these  points,  so 

N(x) 

Q(x)  = ©> 

i= 1 

where  N(x)  is  the  random  number  of  scatterings  on  the  segment  [0,  x]. 

iii)  The  random  variables  Aj,  A2,  A3, . . . are  independent  of  each  other  and  have 
the  same  distribution  density  q(x). 

iv)  The  random  variables  ©i,  ©2,  ©3, . . . are  independent  of  each  other  and  of 
Ai,  A2,  A3, . . . and  identically  distributed  with  the  density  cr(9). 

v)  The  distribution  a (6)  is  concentrated  in  such  a small  region  near  zero  that  the 
vectors  0;  and  Q(x)  may  be  considered  as  two-dimensional  vectors  and  integration 
with  respect  to  6 may  be  extended  to  the  whole  plane  R2. 

On  these  assumptions,  the  angular  distribution  density  p(6,  x)  of  the  particle 
passed  the  path  x is  given  by  the  sum 

X 

p(0,  x)  s=  6(9)Q(x)  + a(6)  J Q(x  - x')q(x')dx'+ 

0 

X 

o*o{0)  J Q(x  — x)q  * q(x)dx,Jt 
0 


a * a * a 


X 

(0)  J Q(x  — x')q  * q * q(x')dx' + . . . 


Here  * means  convolution  of  distributions 


X 

q * q(x)  = J q( x — x')q(x')dx' , 
0 

a * a{6)  = J a(6  — 6')a(9')d9' , 


(1) 
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and 

OO 

Q(x)  = J q{x')dx  . 

X 

Represent  (1)  in  the  form 

X 

p(0 , x)  = J Q(x  — x')f(9 , x')dx  . (2) 

o 

The  function  f(0 , a:)  is  the  density  of  collisions  and  can  be  expanded  in  the  collision 
number  series 

f(d,  x)  = S(6)S(x)  + a(9)q(x)+ 

+cr  * a (9)q  * q(x)  + a * cr  * &{8)q  * q * q(x)  + . . . = 


= £ff(nWn)(a0-  (3) 

n= 0 

where 


= (5(0)^  cr^^(^)  = (t(6),  <7^(6)  = a * . . .a(6),  n>  1 

' v ' 

n 

and  so  are  q(n\x).  It  is  easy  to  see  that  the  distribution  (3)  is  a solution  of  the 
master  equation 


X 

f(0,x)  = 6(9)5(x)  + J dx'q(x—x')  J d0'a(9  — 9')f(d'  ,x').  (4) 

0 R2 

Formula  (2)  together  with  master  equation  (4)  gives  a complete  description  of  the 
problem. 


3 The  scattering  theory  for  a homogeneous  medium 

The  ordinary  theory  of  multiple  scattering  appears  as  a result  of  the  addition  of 
two  further  postulates  to  those  mentioned  above: 

vi)  The  mean  square  angle  of  a single  scattering  has  a finite  value 

OO 

(02)  = J a{9)62d9  = 2tt  J <r(9)93d9  < oo. 

R2  o 

vii)  The  random  variables  Xi,Xi  + X2,Xi  + X2  + ^3, . . . form  the  Poisson 
ensemble,  i.e.  the  increments  X\,  X2,  X2, . . . are  distributed  according  to  the  expo- 
nential law 


q(x)  — pe  p > 0. 


(5) 
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Inserting  (5)  into  (2)  and  (4)  we  obtain  two  expressions 

X 

p(0 \ x)  — J e^x  f(9 \ x')dx\ 

o 

X 

f(0,  x ) = 5(9)6( x)  + J d0’a{9-  O')^*  J e^'  f {O' , x')dx' 

R2  0 

which  after  simple  transformations  take  a form  of  the  integro-differential  kinetic 
equation  (the  Boltzmann  equation) 

= -fJ.p{0,  x)  + p J d,9'cr(d  - d')p(0',x) 

R2 

with  the  initial  condition 

p(0,G)  = 6(0). 

In  the  large  a;-asymptotics,  the  equation  reduces  to  the  ordinary  diffusion  equation 

dp(0,x) 

~dx  ~ = D\7gp(0,x)  (6) 

with  D = p(@'2)/2.  It  is  supposed  here  that  a(0)  is  an  axially  symmetric  function 
depending  on  |0|  = 0. 

4 The  scattering  theory  for  a fractal  medium 

A passage  from  a homogeneous  medium  to  a fractal  one  is  performed  by  the  re- 
placement of  postulates  vi)  and  vii)  with  the  following: 

vi) *  The  angular  distribution  of  particles  scattered  on  a single  fractal  cluster  is 
characterized  by  the  power  law  1 

J a(9')d.0'  ~ A9~a,  0 ->  oo,  0 < a < 2.  (7) 

\§'\>e 

vii) *  The  free  path  distribution  in  a fractal  path  medium  has  a long  power  tail 
too  2: 

Q(x ) ~ Bx~P , x — > oo,  0 < /3  < 1.  (8) 

Note  that  the  mean  square  angle  in  the  first  case  and  mean  free  path  in  the  second 
one  are  infinite. 

The  information  given  by  (7)  and  (8)  is  enough  to  find  the  asymptotic  behaviour 
of  the  distribution  p(9,  x)  at  large  distances  x.  To  do  this  we  apply  the  Fourier- 
Laplace  transform 

OO 

p(k,  A)  = j d9  J dxeiB-Xxp(0,x), 

R 0 
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where  A;  is  a two-dimensional  vector  and  k6  is  the  scalar  product.  Formulas  (2)  and 
(4)  take  the  following  form: 

p{k,\)  = [1  -q(X)]f(k,  A)/A, 


f(k,  A)  = 1 + q(X)a(k)f(k , A). 


Combining  them  we  get 


p(k,  A) 


1 - g(A) 

A[1  - g(A)o-(£)] 


This  is  the  Montroll- Weiss  result  for  random  walk  with  trapping  in  a two- 
dimensional  space  6 (see,  also  7 and  8). 

According  to  the  Tauberian  theorem  9,  the  relations  (7)  and  (8)  determine 
behaviour  of  transforms  a(k)  and  q( A)  at  small  arguments: 


1 — cr(k)  ~ aka , k — » 0, 


a = 


2~°T(1  — a/2) 
r(l  + a/2)  ’ 


1 — q(\)  ~ b\P , A — > 0,  6 = r(l-/3)B. 

In  their  turn,  these  expressions  make  it  possible  to  find  the  asymptotic  form  for 
p{k,  A): 

p(k,  A)  ~ pas(k,  A)  = ^ + cka , A — > 0,  fc  — » 0 (9) 

with  C = a/b,  whence 


Pas  ($i ; 


(2' 


IdxI 

L R2 


dk 


\P~ 


\P  + Cka 


Xx  


= (cafiy^v^ieicx*)-1/*)  (io) 

where 

J dx  I <u> 

L'  R2 

On  the  other  hand,  Eq.  (9)  rewritten  in  the  form 

A tposik.  A)  = -Ckapas{k,  A)  + \P-\ 

can  be  considered  as  the  Fourier-Laplace  transform  of  the  fractional  differential 
equation.  In  terms  used  in  the  book  10  it  has  the  form: 

DPxO+pas0,  x)  = - C(-Ae)a/2pas(e. , x)  + r^_p}S(x)  (12) 

with  the  initial  condition 

Pas(6,  0)  = 5(9). 

When  ot  — * 2 and  /3  — > 1 the  equation  reduces  to  the  ordinary  diffusion  equation 
(6)  with  D = C. 
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5 Results 


The  integral  (11)  giving  the  solution  of  the  equation  (12)  via  (10)  may  be  repre- 
sented in  two  forms  11 . The  first  of  them  is 

CO 

'l4“l/3)(0)  = J g^a\0xp/a)g^i(x)x2p/adx  (13) 

0 


where  gf'*  ( x ) and  (0)  are  the  one-dimensional  one-sided  stable  density  and  the 

two-dimensional  axially  symmetrical  stable  density  determined  by  their  Laplace 
and  Fourier  transforms  respectively: 


J e~Xx  gf\x)dx  = 

0 

and 

J eimg{2a)(0)d0  = e~ka. 

R2 

Stable  distributions  are  described  in  detail  in  11 ,12. 

The  second  form  uses  the  //-functions  - generalized  hypergeometrical  functions, 
also  called  the  Fox  functions  13 , 

irpq  (al’  ai)  ■ • ■ (*p>  ap) 

mn\  (h,Pi)-.-(bq,Pq) 


where  m,  n,p,  and  q are  integer  numbers  such  that  0 < n < p and  1 < m < q\  a, 
and  f3j  are  positive, 

sjk  = (bj  + k)/(3j 


Cjk 


n r(bi  - Pisjk)  n r(!  - ai + aisjk) 

i=i, w /= i 

n r(l  - 6,  + f3lSjk)  n r(a;  — onSjk) 

l=m+ 1 /=n+ 1 


In  terms  of  //-function,  the  angular  function  F),  'p  {0)  looks  as  follows 


F 


«-!(! 


2+a 


H< 


12 

32 


(-1,  l/a)(-a/2,  /?/2)(l  - a/2, 1/2)  \ 
(0,l/a)(— 1,1/a)  ) 


a < 


*ia'p)(0)  = 


1 


r H. 


21 


?r(a0)2  23  v 2 

When  a = 2 and  (3=1 


(l,l/a)(l,/?/a) 
(l,l/a)(l,l/2)(l,l/2  ))' 


1 < a < 2 


Sp^x)  = 5(x  - 1) 
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5f)(^)  = ^exp{-02/4)  (14) 

and  the  density  (13)  reduces  to  the  ordinary  result  for  a homogeneous  medium 
pas{0,x)=  -^-j-exp{-62/(4Cx)} 

with  the  diffusivity  C.  This  is  the  two-dimensional  Gauss  distribution  with  the 
mean  square  angle 

(02>  = J Vasil  x)d2de  = 4 Cx,  (15) 

R2 

describing  a normal  diffusion. 

When  a < 2 and  (3  = 1 we  obtain 

Pasilx)  = (Cx)-2/ag%*}  (9(Cx)-1/a). 

The  mean  square  angle  is  infinite  now  (due  to  the  property  of  the  density  g^), 
but  we  may  take 

69  = {Cxf/a 

as  the  measure  of  the  angular  distribution  width.  This  formula  reveals  a more 
rapid  broadening  of  the  angular  distribution  than  in  the  normal  case  (15).  We 
see  here  a superdiffusion  regime  considered  in  1 (a  more  detailed  consideration  for 
TV-dimensional  space  is  performed  in  14). 

When  we  have  the  opposite  case,  a = 2 and  f3  < 1,  the  mean  square  root  of  the 
scattering  angle  9 exists  and  is 

^ = <> <L 
This  is  a subdiffusion  15 . 

In  the  general  case  the  rate  of  broadening  the  angular  distribution  is  determined 
by  the  exponent  P/or. 

69  = C1/axp/a. 

Thus  the  multiple  scattering  reveals  the  subdiffusion  asymptotic  behaviour  if  P/a  < 
1/2  and  superdiffusion  one  if  P/a  > 1/2.  Both  the  regimes  are  covered  by  term 
’’anomalous  diffusion”  (see  8,16  and  references  therein). 

However,  not  only  the  rate  of  broadening  the  angular  distribution  changes  by 
passing  from  homogeneous  medium  to  fractal  one:  the  very  shape  of  the  density 
changes  too.  Whereas  ’ (@)  is  simply  the  two-dimensional  Gaussian  distribu- 
tion, the  densities  ^(0),  a < 2,  are  two-dimensional  symmetrical  stable  densi- 
ties with  characteristic  exponent  a.  They  have  a dome-shaped  form  near  zero  with 
the  maximum  value 

\k(,Q,1)(0)  = r(1  + 2/a), 

47T 
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long  tails  of  a power  type 


(*) 


T(1  + a/2) 


2n 


\(a,K/2)9 


-a-2 


9 — » oo 


and  no  moments  of  the  order  > a (see  10).  On  the  contrary,  the  densities  ^>^’^(6) 
(3  < 1,  diverge  logarithmically  at  zero 

,(2,/3)/^\  . Iln(9l 


27rr(l  - /?)  ’ 


e^o. 


rapidly  falling  tails 


«(3/3/2  — 1)/(2— /3) 


x(0/2)-2O-W(2-«exp{-(2-^/(2^)(0/2)2/(2-«},  0 - oo  (16) 

and  finite  moments  of  all  positive  orders 

J 92n^\e)de  = 4-[y^+1i))]2. 

R2 

Note  that  (16)  reduces  to  (14)  when  /3  — > 1. 

In  the  case  a = 1,  /3  = 1/2,  the  density  M>2 (^)  is  expressed  in  terms  of  the 
incomplete  gamma  function: 

^•1/2\e)  = (47r)-3/2e02/4r(-l/2,6>2/4). 
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